Dark Matter Halos with Cores from Hierarchical Structure Formation 
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We show that dark matter emerging from late decays [z < 1000) produces a Unear power spec- 
trum identical to that of Cold Dark Matter (CDM) on all observationally relavant scales (> 0.1 
Mpc), and simultaneously generates observable constant-density cores in small dark matter halos. 
We refer to this class of models as meta-Cold Dark Matter (mCDM), because it is born with non- 
relativistic velocities from the decays of cold thermal relics. The constant-density cores are a result 
of the low phase-space density of mCDM at birth. Warm dark matter cannot produce similar size 
phase-space limited cores without saturating the Lya power spectrum bounds. Dark matter domi- 
nated galaxy rotation curves and stellar velocity dispersion profiles may provide the best means to 
discriminate between mCDM and CDM. mCDM candidates are motivated by the particle spectrum 
of supersymmetric and extra dimensional extensions to the standard model of particle physics. 

PACS numbers: 95.30.Cq, 95.35.-|-d, 98.80.-k 
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INTRODUCTION 

In the standard cosmological model, the bulk of the 
matter in the universe is in the form of Cold Dark Mat- 
ter (CDM). Structure in CDM is built up hierarchically, 
from initial density perturbations that are adiabatic and 
nearly scale- invariant [l]. CDM provides a compelling 
description of the universe on large scales, but it has not 
conclusively passed all observational tests on small scales. 
Two vifell-known shortcomings are that simulations pre- 
dict Milky Way-sized galaxies should have nearly an or- 
der of magnitude more dark matter substructures than 
observed satellites |2|, and density profiles of low-mass 
galaxies tend to be less cuspy than those seen in simula- 
tions [1, 0, Ml @1 and are often well-fit by constant-density 
cores. While the former can be alleviated both by astro- 
physical and cosmological [1, @] methods, the latter, if 
ultimately true, may require a more drastic modification 
of the CDM paradigm 10]. 

Warm dark matter (WDM) provides an alternative 
that may alleviate the problems on small scales 11, T3 |. 
Indeed, we expect constant density cores in the centers 
of WDM halos because the initial phase-space maximum 
can never be exceeded [l3|. An additional consequence 
of the large WDM particle velocities is that its free- 
streaming length is large. This suppresses the power 
spectrum on small scales. One problem with WDM as 
an alternative to CDM, however, is that models with ve- 
locities large enough to produce observationally-relevant 
cores are in conflict with measurement of the Lyman-a 
forest power spectrum 1J| . Indeed, Lyman-a forest anal- 



yses suggest that the linear density fluctuation spectrum 
is nearly identical to the CDM prediction on scales larger 
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in the halos of low-mass galaxies. In this paper we in- 
vestigate such a dark matter model, which we refer to as 
"me to-Cold Dark Matter" (toCDM). We define mCDM 
as particles that emerge relatively late in cosmic time 
(z < 1000) and are born non-relativistic from the decays 
of cold particles. As we will show, the mCDM model pro- 
vides both large phase-space cores and CDM-like power 
spectra on > 0.1 Mpc scales, consistent with the most 
stringent measurements of the Lya forest power spec- 
trum [3 • Specifically, mCDM is born late enough 
that the free-streaming scale is small, but the velocity 
dispersion is large enough to give rise to a reduced phase- 
space density. Even though mCDM velocities are non- 
relativistic, the redshifted velocities today are larger than 
the corresponding velocities for WDM. Previous studies 
of similar models considered decays with shorter lifetimes 
~ M^,/M„?.„„t. (of order a month) to relativistic daugh- 
ter particles [20|, |21|. In these models the relationship 
between the phase-space density and cut-off scale in the 
power spectrum is similar to that of WDM. 



It is thus interesting to consider if any modification to 
CDM can reproduce the pattern of hierarchical structure 
formation, and also create large phase-space limited cores 



Dark matter in the mCDM class arises in many ex- 
tensions to the standard model of particle physics. For 
example in supersymmetric models we can consider a 
sneutrino decaying into a neutrino and gravitino with a 
lifetime r ~ 3.6 x 10^ s (100 GeV/Am)4m(5/TeV 
where Am is the mass difference between the sneutrino 
and gravitino, and ttIq is the gravitino mass. A simi- 
lar relation holds for Kaluza-Klein WIMPs decaying to 
gravitons in theories with extra dimensions. For exam- 
ple, we can consider a 10 TeV gravitino produced from 
sneutrino decays with a lifetime t ~ 5 x 10^^ s, where 
the velocity imparted to the gravitino is vjc ^ 10~^. 
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DECAY KINEMATICS 

We assume that the decays take the form CDM 
toCDM + £, where the decaying CDM is non-relativistic, 
toCDM is the dark matter today, and £ is a hght neutral 
particle. After the epoch of decay, ad, the velocity of 
the daughter particle will be w ~ {pcm./'m){ad/a), where 
771 is the 777CDM mass and Pcm is the center of mass 
momentum of the daughter particles. We define the av- 
erage primordial phase-space density for the 777CDM as 
Qp = p/(T^, where p is the mean density and a is the 
velocity dispersion. Integrating over the phase-space dis- 
tribution function (presented below) we find [20| 
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(1) 



where a = 1.0(0.8) for decays in the radiation (matter)- 
dominated era. Here, and for the rest of the discus- 
sion, the units of Q are M0 pc~^(kms^^)^'^. Note that 
for the typical cases we consider, {mcdm ~ 'ni)/mcdm — 

Pcm/rricdm ^ 10"^. 

In the case of late decays, the power spectrum is gov- 
erned not only by Qp but also by the epoch of decay, ad- 
We can understand this heuristically in terms of standard 
free-streaming arguments. In this approximation, power 
on a comoving scale smaller than a streaming scale, 
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FIG. 1: The ratio of the power spectrum to that of standard 
ACDM for a variety of different models. The red curves show 
A+mCDM models with r = 5 x 10^^ s and Qp = lO"'^ (long- 
dash) and Qp = 10~^ (short-dash). The solid black curve is 
a A+mCDM model with r = lO" s and Qp = 10"*^. Also 
shown are limits on the warm dark matter power spectrum 
from the analyses of [3l (dotted, right) and [T^] (dotted, left). 
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will be erased. Here r = t{ad) is the decay lifetime, to 
is the current epoch, and we have assumed that the de- 
cay is non-relativistic. If the decay were to occur well 
before matter-radiation equality, ad < a^q, we would 
have Xps — Pcmad/m\ii{aeq/ ad)ail^Q.hH^^ , and the 
filtering would be driven primarily by the phase-space 
density variable PcmO.d/'Tn, and only logarithmically on 
the decay lifetime through ad- If we consider decays 
in the matter-dominated regime, then the power spec- 
trum filtering scale depends on both the decay epoch and 
the phase-space density: \ps — {jpcmO'd/'m)a^^^^Q-bHQ^ 
(~ 0.05 h~^Mpc). Thus for long lifetimes (a^ > Ogg), the 
free-streaming scale is suppressed even if Qp is relatively 
small. For the models we consider, the universe does not 
re-enter a radiation dominated phase due to the decay. 



POWER SPECTRUM 

In order to track the evolution of fluctuations we will 
need the ttiCDM phase-space distribution. Following [2^, 
24 | , we find that 
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The variables aq = q/pcrm and tq = t(aq) are charac- 
teristic epochs and times for each value of the comoving 
momentum q- In the limit of small Pcm/Tn, we can write 
all of the perturbations to the phase-space distribution 
in terms of the density and velocity perturbations. For 
the decaying CDM, the evolution of these quantities is 
identical to the standard non-decaying case [2^. This 
results from the fact that the decays act to remove the 
same number of particles from every region of space. 

To compute the power spectrum we use the full hi- 
erarchy of equations governing the perturbation of the 
phase-space distribution. We use the concordance cos- 
mology [2^, and assume that the universe is made flat 
by a cosmological constant (A), and that the dark mat- 
ter is in the form of 77iCDM. We can gain insight into 
the effect of the decays by applying energy conserva- 
tion and considering a simpler set of equations in the 
limits where the 777,CDM equation of state and sound 
speed is negligible. For large scales, fc ^ 0, we have 

y' = 5" - = -2//t(1 - 277igpcdm/TOcdmP)' ^hcrc 5 is 
the 771CDM density perturbation and the derivatives here 
are with respect to cosmic time, related to the conformal 
time, rj, by drj — dt/a. We have defined 777q = rn^^^ — m?. 
For small a, y is driven to zero as and is exponen- 
tially damped for a > ad- Thus the ttiCDM perturbation 
'catches up' with the CDM perturbation almost immedi- 
ately, and the power spectrum retains the form of stan- 
dard ACDM on scales larger than the horizon at decay. 
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We expect departures from a ACDM power spectrum 
approximately when k^c^ri^ ^ 1, where is the mCDM 
sound speed. For decays in the radiation-dominated 
regime, this damping scale is given by PcmOd/w, and 
for decays in the matter-dominated regime the scaling 
is (j)cmO'd/'>Ti)T^^^'^ ■ Thus the damping scale for matter- 
dominated decays is set by Qp and r, while for decays 
during the radiation era the scale depends only on Qp. 
As one would expect, these are the same scalings we ob- 
tained earlier with our simple free-streaming analysis. 

In Figure [1] we show the ratio of the power spectrum 
for a variety of A-^mCDM models to ACDM. The dotted 
lines are warm dark matter models that correspond to 
the bounds from recent analyses of the Lyman-a forest 
[illii. The model with r = 5 x lO^^ s and Qp = 10"^ is 
potentially inconsistent, though these bounds can clearly 
be evaded with a Qp = 10~^ model. An example model 
with a longer lifetime of 10^'* s and Qp — 10~^, however, 
is clearly consistent with both bounds, which is what we 
predict from the above scalings. Further detailed model- 
ing of the Lyman-a power spectrum on small scales thus 
provides an excellent probe of viable mCDM models. 



PHASE-SPACE CORES 

We can compare the primordial phase-space values Qp 
to the phase-space densities deduced from measurements 
of galaxy rotation curves. Among the highest-resolution 
rotation curves for low-mass spirals are those presented 
by [3] who use two-dimensional Ha and CO velocity 
fields. When the dark matter halo components of these 
galaxies are fit to cored density distributions, the aver- 
age central phase-space densities are 2 x 10-6 1,113. 
cannot be smaller than this value. 

In order to connect Qp to the maximum phase space 
density in a halo, Qq, we must consider halo formation. 
We apply a general argument for coUisionless collapse 
that is much stronger than the Tremaine-Gunn bound 
for distributions with maximum fine-grained phase space 
density much larger Qp. Specifically, the Excess Mass 
Function (EMF), defined in 1271 . must always decrease. 
Following the approach of [2(J|, we calculate the EMF 
of primordial mCDM using Eq. [3] and demand that its 
value for all phase-space densities be greater than the 
EMF of the dark matter halo profile. This procedure 
gives us the maximum central phase space density of the 
dark matter halo profile that is consistent with a given 
primordial momentum distribution. 

The procedure depends on the details of the decay pro- 
cess as well as the assumed dark matter halo profile. We 
will restrict our attention here to the Milky Way satel- 
lite galaxies, and this motivates_a general parameteriza- 
tion for the dark matter halo 



that interpolates be- 
tween a truncated NFW profile and a profile with a core 
- p(r) ~ exp(— r/6)/(c -\- r)/{r + a)^. For parameters 



that are broadly consistent with the Milky Way dwarf- 
satellite galaxies we find numerically that Qo ^ 30(5p for 
Qp = 10-5 and go < 50Qp for Qp = 10-^. 

The above procedure also provides an estimate of the 
minimum core size. We find that the minimum value of 
c ranges between 300 to 800 pc for Qp = 10^® depending 
on the value of a (restricted to be greater than c) for a 
halo with mass within 1 kpc of IO^Mq (which results in 
a central stellar velocity dispersion ~ 10 km/s for a typi- 
cal satellite dwarf galaxy) . The lifetime is consistent with 
the model in Figured] A decrease in halo mass to IO^Mq 
increases minimum c range to 1.5-2.5 kpc, while increas- 
ing mass to lO^M© decreases it to 150-300 pc range. The 
minimum c scales approximately as l/^/Qp in this range 
of parameter space. 

We caution that the above exercise only provides us 
with the minimum possible core size. The merging pro- 
cess is expected to increase this core size [2I] and detailed 
simulations are required to make further predictions. In 
the above example, we have allowed the inner core (c) 
and the overall scale radius of the halo (a) to vary freely 
of each other. However, mergers will certainly correlate 
the two scales. Also, since halos get built out of smaller 
halos, the excess mass function analysis should be applied 
to the smaller building blocks before merger. 

Hierarchically-formed dark matter halos from CDM 
particles without an appreciable phase-space maximum 
are known to have power-law phase-space profiles, 
Q(r) = p{r)/a^{r) = Q_2(r/r„2)~", with a ~ 1.875 
[29|. Here we have normalized the Q profile at the ra- 
dius where the log-slope in the density profile is —2. We 
can use this result also to estimate the size of the core 
radius, Tcorc, that would be imposed for a given central 
phase-space density. A minimum estimate for the core 
size will assume that the power-law continues uninter- 
rupted until the central Q is reached, Q{r = rcoic) = 
Qq. We must first determine the normalization con- 
stant Q-2- We do so assuming an NFW profile (s^l 
and find Q(r) ~ lO-^^M^^\c„r/l5)-°-^^Hr/R^)-^-'^^^, 
where Mn = (A'/v/lO^^/j-^M©) characterizes the halo 

1 /3 

virial mass, i?v — 134 kpc Mj^{ is the halo virial 
radius, and Cmr = Rv/r^2 — 15 is appropriate for 
low-mass spiral galaxy halos The quoted power- 

law dependence is approximate but serves to highlight 
that the dependence is expected to be weak. If we 
use Q{r coio) = Qo to derive a lower limit on the core 
size and ignore weak Cyir dependence we obtain rcore ^ 
2.4 X 10-3i?,Mfi°-53g-o.53^ ^ Qo/(10-'5) 

characterizes the Qo dependence. Small spiral galaxy ha- 
los (Mil — 1) with Qo set by the Lyman-a forest limit for 
WDM (Q_6 — 8 X 10'^) can have only small cores rcorc ^ 1 
pc. A toCDM candidate can avoid the Lyman-a forest 
bound with Q_6 ^ 10 and produce larger, dynamically- 
important cores, rcorc ^ 100 pc. As we have emphasized, 
this estimate represents a conservative lower limit on the 
full extent of the core region because we have assumed 
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a sharp break in the phase-space profile. The WDM N- 
body simulation results of [S^l suggest that the cores are 
^ 3 times larger than our estimate would give. 

Note also that the core size is expected to take up a 
larger fraction of the halo as we consider smaller sys- 
tems. Large, soft cores will render these dwarf-size halos 
uite prone to disruption upon accretion into larger halos 
In mCDM we thus expect the predicted substructure 
count to be reduced relative to CDM predictions. 



DISCUSSION 

For the models considered, the net effect of the injected 
relativistic energy during decay may be phrased in terms 
of the effective number of light neutrinos, AiV,y ~ 1.8 x 
10~^-\/r/yrs(pcm/TO). For the = 10 TeV, t — 5 x 
10^^ s model described above, we have AN^ ~ 0.01 which 
would be hard to detect. 

Neutrinos and photons produced along with mCDM 
will arrive unscattered in the diffuse radiation back- 
grounds today. The strongest constraint comes from the 
neutrinos produced directly in two-body decays. The rel- 
evant low-energy neutrino flux limit is < 1.2 cm~^ s 



from Super-Kamiokande for energies 18 — 82 MeV [33 1. 
The above sneutrino model is consistent with this limit 
and in the sensitivity range of future experiments with 
reduced energy thresholds [3J] . Photons may also be pro- 
duced in three-body decays and subsequent hadroniza- 
tion into neutral pions. Using a very conservative branch- 
ing fraction of 10^"^ into photons [HH], we find that the 
photon flux is consistent with observations [36| . 

In conclusion, we have shown if dark matter is pro- 
duced from the late decay of cold relics, it will give 
rise to large cores in low-mass galaxies and alleviate the 
dwarf satellite problem, without destroying agreement 
with the observed Lyman-a forest power spectrum. This 
differs from standard WDM models, which, in order to 
remain consistent with constraints from the Lyman-a 
power spectrum, cannot produce sizeable cores in small 
galaxies. Dynamical studies of nearby galaxies may pro- 
vide the best means to test the mCDM scenario and com- 
pare it to CDM predictions. 
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